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periodic modulating field. The dependence of the revival time on various 
parameters of the driven system is shown analytically As an example of 
application of our approach, we compare the analytically obtained values of 
the revival time for various modulation strengths with the numerically 
computed ones in the case of a driven gravitational cavity. We show that 
they are in very good agreement. 
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In nature interference phenomena lead to recurrences [1, 2]. In quantum 
mechanical evolution, for instance, interference plays an important role and 
manifests itself in quantum recurrences [3]. A quantum wave-packet spreads 
all over the available space after a few classical periods following wave me- 
chanics and collapses. However, due to quantum dynamics it rebuilds itself 
after a certain evolution time which is larger than the classical period. In 
one-dimensional systems this phenomenon is well studied and is known as 
quantum revivals and fractional revivals [4, 5, 6, 7, 8, 9]. In higher dimen- 
sional systems, which exhibit classical and quantum chaos, this phenomenon 
has been numerically observed earlier [10, 11, 12, 13, 14]. In particular, the 
quantum revivals occurring in driven gravitational cavities [15, 16] have been 
investigated analytically and numerically. 

In this Letter we provide a general analytical prescription to the revival 
phenomenon occurring in any periodically-driven time-dependent systems. 
By using semi-classical secular theory [17, 18, 19], we derive a simple relation 
which allows us to calculate the quantum revival time in the presence of 
the external periodic modulation. For the sake of concreteness, we apply 
our general result to the dynamics of atoms in a modulated gravitational 
cavity [20] which is accessible to laboratory experiments [21]. Hence we 
argue that dynamical revivals are generic to all periodically-driven explicitly 
time-dependent systems, so far as resonances prevail. 

Let us consider a one-dimensional system driven by an external periodic 
field. In order to calculate its quasi-energy, we consider the nonlinear reso- 
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nances of the time-dependent system [17, 18]. We denote the eigenstates and 
eigenvalues of the corresponding time independent system by \n) and E n , re- 
spectively, so that Ho\n) = E n \n), where H is the unperturbed Hamiltonian. 
We consider the evolution of a wave-packet in an arbitrary one-dimensional 
potential in the presence of an external periodic field. The Hamiltonian of 
the driven system in dimensionless form can be expressed as 



H = H + A V(x) sint, 



(1) 



where A is the dimensionless modulation strength and V(x) defines the cou- 
pling. 

In order to study the quantum nonlinear resonances of the system we 
make the ansatz [17, 18] that the solution of the Schrodinger equation cor- 
responding to the Hamiltonian (1) can be written in the form 



\m) = J2 C n(t)\n} expl-i 



E r + (n-r)^ 



k 



(2) 



where k denotes the effective Planck's constant, and E r is the average energy 
of the wave-packet in the Nth. resonance. 

We substitute Eq. (2) into the Schrodinger equation and project the result 
onto the state vector (m|. This leads to 



ikCn 



k 

E m -E r + (m-r) — 



X 



e i(n~m+N)t/N _ i(n-m-N)t/N 



C n (t) , 



(3) 



where V n>m = (n\V(x) \m) are the matrix elements oiV{x). In Eq. (3) we drop 
the fast oscillating terms and keep only the resonant ones, that is, n = m±N. 



Moreover, for large m we may take [17, 18] V m , m+ N ~ V m ,m-N = V. Hence, 
Eq. (3) reduces to 

k 



ik C m — 



E m - E r + (m - r) 



C m (t) 



+-prr (C m +N — Ctti-n) ■ (4) 
Li 

In view of the slow dependence of the energy E m on the quantum number 
m around the initially excited level r in a nonlinear resonance, we expand 
the energy E m up to second order. Hence, the equation of motion for the 
probability amplitude C m takes the form 

mC m = (m-r)^E' r -^jC m 

+ \{m - r) 2 E';C m + ^ (C m+N - C m _ N ) . (5) 

For the exact resonance case E' r = k/N, and as a result the first term on the 
right hand side vanishes. We may write this equation in angle representation 
by introducing the Fourier representation of C m as 

1 2w 

C m = —jgi^e-^-^dcp 
o 

2N-n 

= 2A% / 9(0)e- l{m - r)9/N d9 . (6) 
o 

This particular choice of the Fourier representation of C m yields a Schrodinger- 
like equation, ikg(9,t) = Hg(9,t), where the Hamiltonian H(9) reads 

N 2 E" d 2 

n{9) = — ^w + XVsin{eh (7) 
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and is independent of time. To obtain this equation, we have assumed that 
the function g(9,t) has a 2iV7r-periodicity in ^-coordinate. 

Due to the time-independent behavior of 7~C(9), we can write the time 
evolution of g as g(9,t) = g(9)e~ i£t ^ . On substituting this representation in 
the equation of motion for g{9, t), we are left with the rather simple equation 



N 2 E'J d 2 



r 



-£ + XVsm(9) 



~g{9) = . (8) 



2 d9 2 

With the change of variable 9 = 2z + tt/2, we can reduce it to the standard 
Mathieu equation 

" d 2 



t „ + a — 2qcos(2z) 



~g{z) = . (9) 



Here, a = 8£/N 2 E", q = 4\V/N 2 E", and g((p) is a 7r-periodic function. 

Through this procedure we have mapped our initially unsolvable time- 
dependent Schrodinger equation onto the known Mathieu equation [17, 18] 
under the assumption of exact resonance. The Floquet solution of the Math- 
ieu equation [19] can be written as g v {z) = e lzu P u (z), where P u are the 
Matthieu functions [22]. Since we require a 7r-periodicity of g, we are forced 
to take only the even values of the index, i.e., v — 2k/ N. 

Hence, the Floquet quasi-energy eigenvectors [19] can be written as 

W))=e^Mt)), (10) 

where and \uk(t)) are defined as 

N 2 E" 

£ k = — ^(q), (11) 



27T 

\ Uk (t)) = J_y^ e int/N ( dve ivN^2 e -r(n-r)v/2 

27r n 

xPu(k) \n) . (12) 

In this way, we have obtained an approximate solution for a nonlinear reso- 
nance of our explicitly time-dependent system. 

In order to check the correctness of our result, we study the case of zero 
modulation strength, that is A = 0, for which q — 0. The corresponding 
value for the Mathieu characteristic parameter becomes a v t n \(q = 0) = A(n — 
r) 2 /N 2 . This reduces the quasi-energy given in Eq. (11) to 

£n(q = 0) = ^-(n-r)\ (13) 

which is indeed the energy of the un-modulated system expanded around 
n = r for the exact resonance case. This is obtained by considering k = n — r, 
and therefore we may express u n (q = 0) = 2(n — r)/N. 

The eigenfunctions given in Eq. (10) form a complete set of basis vectors. 
Therefore, we can write the propagated wave-packet at any later time in 
terms of the quasi-energy eigenstates, such that 

WM)}=E£ne-^kW>, (14) 

n 

where £ n describes the probability amplitude in the nth state. 

In order to investigate the dynamical revivals, we calculate the auto- 
correlation function between the initially excited wave-packet and the wave- 
packet \ip(z,t)) after a certain evolution time t, viz. 

(m\m) = T,\u 2 e- i£ - t/Js ■ (is) 

n 
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We take into account that the wave-packet is initially centered around n = r 
in the energy representation. Therefore, by expanding the quasi-energy of 
the system around the resonant level r, we can write 

(m\m) = E IU 2 exp {-i p°) + (n - r)w« 

n 

+ (n-r)V 2 ) + ---]t} , (16) 

where u;^ = (j\kY 1 £^\ n=r and £^| n=r denotes the jth derivative of £ n 
calculated at n — r. 

We then compare the auto-correlation function for the driven system, as 
given in Eq. (16), with the auto-correlation function for the corresponding 
one-dimensional unperturbed system [5, 15]. This comparison helps us to 
identify the classical period for the driven system as 

T - 2n 2nk (17) 

On | n =r 

and the quantum revival time as 

T ^ -J2) = ■ ( 18 ) 

^ On \ n =r 

Let us now consider the case of small perturbations, i.e. A<1. We may 
therefore use the expansion for the Mathieu characteristic parameter a v [22] 
only up to second order in q. Then Eq. (11) becomes 

Sn = T ( z/ ' + T^i) +0(g4) ' (19) 

and a simple expression for the revival time in the driven potential can be 
derived from Eqs. (18) and (19). 
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In general, the initially excited wave-packet is away from the center of 
the resonance and the condition E' r = k/N is not satisfied. For the sake 
of clarity, in the above discussion we have only considered the case of the 
exact resonance. However, taking into account the non-resonant situation, 
the revival time for a time-dependent system can be calculated in the general 
case and reads [15] 



T\ = Tn 



2\Ej \N 



3i£ 



(20) 



Here, T = Airk / E" denotes the revival time of the initial wave-packet in the 
un-driven potential and the index v n is defined as 



n N [ )+ NE'J 



(21) 



where the second term on the right hand side is arising due to the non- 
resonant case. Substituting the value of v n from Eq. (21) calculated at n = r 
in Eq. (20), we obtain the general expression for the quantum revival time 
in a periodically driven system as 



T X = T 



1 - 



1 fXV\ 3fi 2 + N 2 /4 



(22) 



2\E»J (fi 2 -N 2 /A) 3 

where \x = ^ Er ~J/ N ^ . Equation (22) constitutes our main result: it expresses 
the modification of the quantum revival time in the presence of an external 
periodic driving field as a function of the modulation strength A and of the 
other characteristic parameters of the system. 
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The present approach is valid for small modulation strengths. In the 
weak binding potentials, for which the level spacing decreases as we increase 
the quantum number, the nonlinear resonances disappear after small mod- 
ulations, as observed numerically for gravitational cavity [23, 24] and also 
for hydrogen atoms. As a consequence, this technique is reasonably good 
for these kind of potentials. For tight binding potentials, for which the level 
spacing increases as we go up in the energy, we find nonlinear resonances 
for higher modulation as well. In this kind of potential we need to consider 
higher-order terms in the expansion of the Mathieu characteristic parameter 
in order to correctly calculate the quantum recurrence time. 

As an application of our method, we consider the dynamics of atoms in a 
gravitational cavity [20] in the presence of an external periodic field [25] . In 
this case the atoms move in a potential x + Voe~ KX , where the linear term is 
due to the gravitational potential and the exponential term comes from the 
evanescent-wave field which can be obtained by total internal reflection of a 
laser light field on the surface of a glass prism [20, 21, 27, 28, 29, 30]. In the 
approximation of a triangular well potential, the coupling constant for the 
iVth resonance is given by 

2E 

V={m\x\m±N) = - - Q(m _ 2)]2 , (23) 

which in the limit of large m becomes V = —2E N /N 2 ir 2 . Substituting this 
expression into Eq. (22), it is possible to calculate the time of revival for a 
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wave packet initially excited with average energy E r around the iVth reso- 
nance, which reads 



where a = (E N /E r ) 1 ^ 2 and a = a 2 k/AE r . If the initial energy is large 
enough, that is E r ^> 1, we may neglect a 2 with respect to (1 — a) 2 in 
Eq. (24), obtaining the simpler formula 



This expression allows us to compare the analytically calculated revival time 
with the numerically computed ones. 

We may prepare a wave packet by trapping and cooling cesium atoms 
in a MOT [20, 21, 28]. In laboratory experiments the driving frequency 
can be changed from to 2MHz, see Ref. [21]. Hence if we select a fre- 
quency of uj = 2nx0.93KHz out of this domain, for cesium atoms of mass 
M = 2.2xl0~ 25 Kg, the dimensionless Plank's constant is AT = 1. Moreover, 
by using effective Rabi frequency Q e ff = 2rcx3.72KHz and the steepness 
0.57/wn, we have V = 1 and k — 1. In our numerical calculations we use 
these values of the dimensionless parameters. In order to measure the quan- 
tum recurrence time, we place the atomic wave packet at z = 29.8/im for 
which E r = 104.1, and z = 20.1/im, for which E r = 70.28 above the atomic 
mirror. 




(24) 




(25) 
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Our numerical results show a complete qualitative and quantitative agree- 
ment with the analytical results. In fact, from Eq. (25) we learn that (i) the 
revival time changes quadratically as a function of the strength A of the exter- 
nal modulation field, and (ii) the revival time depends inversely on the square 
of the initial average energy of the wave-packet. Our numerical investigation 
confirms both these analytical results. 

In Fig. 1 we compare the quantum revival times calculated from Eq. (25) 
(dashed lines) with those obtained numerically (solid lines) for an atomic 
wave-packet bouncing in a modulated gravitational cavity. We see from the 
numerical data that the revival time displays a quadratic dependence on the 
strength A of the external modulation. Moreover, the change in the revival 
time is just about 8% for A ranging from to 0.25 when the initial average 
energy E r = 104.1, (solid line with circles), whereas it increases to almost 
40% for the same range of A when the initial average energy E r = 70.28 
(solid line with squares), in agreement with the inverse square dependence 
of T\ on E r . 

In summary, we have presented a general approach to the investigation 
of the phenomenon of quantum revivals in periodically driven systems. We 
have derived a simple relation which provides the quantum revival time as 
a function of the modulation strength. We have finally applied our theory 
to the case of the revival time arising in the dynamics of atoms bouncing 
in a modulated gravitational cavity, obtaining an excellent agreement with 
the exact numerical data. For the set of parameters which we have used in 
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our calculations a modulated gravitational cavity can be realized within the 
framework of presently available technology [15]. Therefore we suggest an 
experiment, designed with the help of these parameters, to test the quantum 
recurrences in the driven gravitational cavity. 

We thank G. Alber, I. Bialynicki-Birula, R. Grimm, W. P. Schleich, and 
F. Steiner for many fruitful discussions. We gratefully thank Prof. Wolfgang 
Schleich for his hospitality at Ulm, where part of the work was carried out. 
FS acknowledges the support from QAU, and MF acknowledges financial 
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